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ABSTRACT 


Conformal mapping was used to obtain the free streamlines and wall pressure coef- 
ficient for jets issuing from two nozzles and striking a surface. The jet nozzle exit 
planes are at an angle to the surface and can be spaced an arbitrary distance apart and 
an arbitrary distance from the surface. The effect of nozzle position and tilt anlge on 
the extent of the flow back between the nozzles is shown. The back flow results from a 
collision between part of the flow from each jet that has been turned by striking the plate. 



TWO DIMENSIONAL INVISCID JET FLOW FROM TWO NOZZLES 
AT AN ANGLE TO A PLANE SURFACE 
by Marvin E. Goldstein and Robert Siegel 
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SUMMARY 

Conformal mapping was used to obtain the free streamlines and wall pressure coef- 
ficient for jets issuing from two nozzles and striking a surface. The jet nozzle exit 
planes are at an angle to the surface and can be spaced an arbitrary distance apart and 
an arbitrary distance from the surface. The effect of nozzle position and tilt angle on 
the extent of the flow back between the nozzles is shown. The back flow results from a 
collision between part of the flow from each jet that has been turned by striking the plate. 


INTRODUCTION 

Single and multiple fluid jets of various cross-sections and in various orientations 
impinging on surfaces are used in a number of fluid mechanical devices and industrial 
heat and mass transfer applications. These include ground effect machines, V/STOL 
aircraft, cooling of steel strip, tempering of glass and drying of paper. Some of the heat 
transfer applications are discussed in reference 1. 

The impingement features of jets are not fully known, especially in the case of 
multiple jets. Some insight can be obtained by considering inviscid slot jets, as these 
can usually be analyzed by employing conformal transformations. An example is refer- 
ence 2 where the flow pattern was determined for two parallel jets striking a plate. The 
potential flow solution is required to compute the boundary layer along the surface needed 
to determine the heat transfer or mass transfer. 

This report presents an analytical solution for the flow field produced by two inviscid 
two-dimensional jets striking a plane surface. The jet nozzle openings are at an angle 
to the surface as shown in sketch (a), and the surrounding fluid is stationary. The study 
was undertaken to provide a qualitative picture of the flow pattern under a vertical take- 
off aircraft with a fan pod in each wing as shown in plan view in sketch (b). Although the 
computed flow is idealized, it is hoped that it will lead to an increased understanding of 





some aspects of the flow field beneath the aircraft. A consideration in the aircraft 
design is the recirculation upward under the fuselage as illustrated in sketch (c). Although 
this can provide a beneficial lifting force it can also provide an undesirable circulation 
around poritions of the aircraft of hot gases carrying debris from the ground (refs. 3 
to 5). If the engines are tilted outward, the recirculation may be reduced, and for a 
sufficiently large tilting there will be no recirculation. In a real jet there is entrain- 
ment along the free streamlines and this will attenuate the velocities of the jet as it 
moves through the surrounding fluid. However, if the jet is only a few nozzle widths 
from the ground, the entrainment will be small and an inviscid solution will provide a 
reasonable guide to the flow distribution. 

Within the assumptions of inviscid slot jets in constant pressure stationary surround- 
ings, the flow patterns obtained here are exact solutions. The basic theory of conformal 
mapping as applied to jets is given in reference 6 among others. The flow configuration 
in the physical plane is found by carrying out an integration of the velocity potential 
multiplied by the reciprocal of the complex conjugate velocity. Conformal mapping is 
used to find the functional relation between the potential and velocity so that the integra- 
tion can be performed. 

The flow pattern will depend on the boundary condition at the nozzle. It is assumed 
here that the flow exits normal to the nozzle opening, which is a reasonable assumption 
for small nozzles or when turning vanes are used. This provides a nonuniform velocity 
distribution leaving the nozzle as a result of the flow interaction with the ground; the 
velocity in the central part of the nozzle is reduced relative to that along the sides of the 
nozzle. Numerical results will be given for the velocity profile at the nozzle exit, the 
free streamlines bounding the flow and the pressure coefficient along the ground. 


SYMBOLS 


A a constant 

A q a constant 

C Q a constant 

Cp pressure coefficient 

I flow region in the £ plane 

J flow region in the W plane 

K complete elliptic integral of the first kind 

K’ defined by K T (k) = K(k T ) 

k modulus of elliptic integral 
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defined as Vl - k 2 
pressure 

ambient pressure outside jets 
functions defined in eqs. (B-4) and (B-5) 
functions defined in eqs. (B-2) and (B-3) 
complex variable in T -plane, £ + i rj 
complex variable in t -plane 
velocity in x direction 
velocity 

velocity at nozzle exit 
average velocity at nozzle exit 
velocity along free streamlines 
velocity in y direction 
complex potential, cp + iip 

coordinates of center of nozzle exit plane nondimensionalized by A 

coordinates in physical plane 

complex variable, x + iy 

tilt angle of nozzle, sketch (d) 

flow region in the T -plane 

width of nozzle 

asymptotic widths of streams flowing to left and to right 

complex conjugate velocity, u - iv 

Theta -function 

coordinates of T -plane 

density 

quantity defined as iK’/2K 
velocity potential 
stream function 
function defined in eq. (3) 
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Subscripts: 

G, H refer to stagnation points, sketch (d) 

r refers to boundary point, sketch (i) 

Superscripts: 

* complex conjugate 

+ variables defined in eq. (B-l) 


ANALYSIS 

The flow configuration as shown in sketch (a) is symmetric about a vertical plane 
centered between the nozzles, hence only half the flow field need be considered. The 
boundaries in the physical plane (z plane) are shown in sketch (d). Since the region outside 
the jets is at constant pressure, the Bernoulli equation shows that the velocity has a con- 
stant magnitude along the free streamlines DE and CF. The direction of the velocity on 
GF and GHE must be along the boundaries, which correspond respectively to the plane of 
symmetry and the ground. The point H is a stagnation point at which the direction of 
the velocity along GHE must reverse. The dashed line is a streamline dividing the 
portions of the flow that go to the right along the ground or to the left into the stream that 
flows back upward. The point G is a stagnation point. Along the nozzle exit CD the 
velocity is specified as being perpendicular to the line CD and hence is of constant direc- 
tion. 


F 



Id) Boundaries in physical plane (z-plane) where z = x + iy. 
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Let f be the complex conjugate velocity u - iv. From the velocity conditions on 
the flow boundaries it is deduced that the flow field in sketch (d) maps into the interior of 
the region I shown in sketch (e) of the hodograph plane (f -plane). (Corresponding points 
in the various planes are denoted by the same letters). The free streamlines map into 
portions of a circle since they have a constant velocity magnitude. The velocity at the 
nozzle exit maps into a cut at constant angle. 

Let W = cp + iip be the complex potential of the flow, where cp is the velocity po-, 
tential and ip is the stream function. The following considerations show that the flow 
field maps into the interior of the region J in sketch (f). Since the velocities are perpen- 
dicular to the nozzle exit, the streamlines (lines of constant ip) must intersect the nozzle 


-iv 



(e) Hodograph plane (£-plane) where £ = u - iv. 
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(0 Complex potential plane (W-plane) where W - <p + Up. 
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exit CD at right angles. Then since the equal potential lines also intersect the stream- 
lines at right angles, we conclude that the nozzle exit CD must be a line of constant cp. 
The stream function ip is constant along the free streamlines DE and CF and along the 
boundary streamlines FGHE. The dividing streamline of the flow must extend from a 
point on the nozzle exit to the stagnation point H. It follows from these conditions that 
the flow being considered maps into the interior of the region J shown in sketch (f) of the 
complex potential plane (W -plane). 

In order to obtain the solution of the problem the known flow configuration in the 
potential plane is mapped into the physical plane by performing the integration 


z 



~dW + constant 
£ 


( 1 ) 


To carry out this integration the functional relation must be known between £ and W. 

To find the relation between £ and W it is sufficient (ref. 6) to map some region 
T of an intermediate T -plane (T = £ + iT?) conformally onto'*' the region I of the hodo- 
graph plane and also conformally onto the region J of the W-plane. If £ and W are the 
functions that will respectively map T onto I and J, it is necessary that £ and W 
map the boundary of T onto the boundaries of I and J respectively. If y is any 
section of the boundary of T such that £(y) is one of the labeled sections of the boundary 
of I shown in sketch (e),W(y) must be the corresponding section of the boundary of J 
shown in sketch (f) . 

A convenient region to work with in the T -plane is a rectangle. As will be shown 
the region I in sketch (e) can be mapped into a rectangle by use of Theta -functions. The 
region J in sketch (f) is a generalized polygon which can be mapped into the upper half of 
the t-plane in sketch (g) by a Schwarz -Christoffel transformation. Then a second applica- 



(g) Intermediate t-plane. 


^The term "onto" includes the region and its boundary. 
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tion of the Schwarz -Christoffel transformation will relate the rectangle in the T -plane 
to the upper half of the t -plane. These transformations will involve elliptic functions. 

Let K = K(k) be the complete elliptic integral of the first kind of modulus k, 

0 < k •< 1. Thus (ref. 7) 



( 2 ) 


and let K' = K'(k) be defined by 


K'(k) = K(k T ) 


where 

k T = ^1 -k 2 

Now let T be the rectangular region of the T -plane shown in sketch (h), with vertices 
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at the points -K + iK’, K + iK’, K + iO and -K + iO. The nozzle exit will be mapped 
onto the bottom side of the rectangle. The free streamlines will be mapped onto the 
vertical sides, and the ground and line of symmetry are mapped onto the top of the 
rectangle. The correspondence can be seen by comparing sketches (d) and (h). 


Mapping Functions 


Mapping function between £ and T planes. - The function £: r — I will now be 
constructed. To accomplish this we will use the holomorphic function Q, studied in 
appendix A (see ref. 6) and defined on T by 


n(T; | r ,k) = 


— (T - |J 

iK'\ 

V 4K r 

2K ) 


Ter 


(T + | - 2K) 

4K r 


(3) 


The quantity i9^(Z|r) is the fundamental Theta -function (See discussion in appendix A). 
We will now demonstrate that the proper mapping function £ is given by 


£(T) = V 0 e 1 K 7r / 2)-Q( ^(T; 4 h , k) y&(T; £ Q , k) Ter (4) 

First equation (A17a) shows that £(T) = 0 only at the stagnation points G and H where 
T = + iK* and ijjj + iK', respectively. Equation (A18) shows that the arguments of 

both f2 functions in equation (4) are 0 along DC since the £2 are positive and real on the 
real axis of the T -plane. Then 


arg £(T) = - - a 

2 


for T e DC 


which represents the nozzle cut in sketch (e). Equations (A20) and (A22) show that 

r T e CF 

U(T)| = V Q for< 

It e ED 


which is the proper behavior along the circular hodograph of the free streamlines in 
sketch (e). Next equation (A24) shows that in traversing the points G and H along the 
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boundary of the rectangle r in a counterclockwise direction (i. e. , along FGHE) the 
argument of £(T) decreases by it/2 at G and v at H, and is constant on the segments 
'FG', GH and HE. This agrees with the hodograph configuration in sketch e. 

Finally, it follows by use of (A24) that arg £(T) = (n/2) - a + 0r/2)[(ljj/K) + l] 

+ (ji /4) [(£ q /K) + l] ; T e GF. On the hodograph of sketch (e) arg f(T) = 3zr/2; T e GF. 
The last two expressions will be compatible if we set 



4K 

7T 



+ 2| h 


(5) 


This combined with the maximum modulus principle (ref. 5) shows that the function 
defined by equation (4) with the condition of equation (5) maps T onto I in such a way 
that the boundary of T is mapped onto the boundary of I. 

Mapping funct ion b etween W an d T planes. - In order to complete the relation 
between £ and W the function W that maps T onto J must be constructed. To 
accomplish this the intermediate t -plane shown in sketch (g) is utilized. An application of 
the Schwarz -Christoffel transformation (ref. 8) shows that the mapping that transforms 
the upper half of the t -plane onto the region J of the W -plane in the manner indicated in 
sketches (f) and (g) is defined by 


™ = iA f 
dt ( 


(t - t H ) 


(1 - k‘ 


! t 2 )Vi - 1 2 


; JU t > o 


(6) 


Now it is shown in reference 9 (p. 359) that the rectangle of sketch (h) is mapped 
onto the upper half t-plane of sketch (g) by 

t = sn(T, k) Ter (7a) 
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where sn is the sine amplitude function. The inverse of this function is defined by 

T = f — t (7b) 

J Q Vl - t 2 Vl - k 2 t 2 


it is common to write sn T in place of sn(T, k) when the modulus k is not important to 
the discussion. 

Using the relations (ref. 9) 


Vl - k 2 t 2 - dn(T, k) 


V i - 1 2 = 


cn(T, k) 


Equation (7b) can be written as 


dt 

dT 


= cn T dn T 


At point H 


t H = sn T w = sn( + iK') = 


H 


H 


k sn £ h 

Equation (6) can now be written in terms of the variable T 

(k sn sn T - 1) 


dW _ dW dt_ =iA 
dT dt dT c 


sn T - sn T 


H 


cn T dn T = iA 


dn T cn T 


dn T 


Ter 


( 8 ) 


Upon integration equation (8) becomes 


W = iA< 


k sn | 


/ sn T dT dT _ f dT _ C ol 
dn T J dn T k’ 


= -iA< 


k sn J 


d(cn T) 
dn 2 T 


+ —f dn(T + K)dT + — 

k' J k’ 
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where the following identities (ref. 9) have been used 


Noting that (ref. 9) 


d(cn T) = -sn TdnTdT 
dn(T + K) = 

dn T 

J' dn T dT = sin~*(sn T) = amT 
dn 2 T = k ,2 +k 2 cn 2 T 


we get 



(9a) 


An alternate form is 

W = - ^ jsn tan -1 ^- cn t) + tan -1 (k' tn T) + I + C Q | (9b) 

where the relation (ref. 7) am(T + K) = tan'^k' tnT) + ( n/2 ) has been used. 

Since the velocity potential is constant across the nozzle exit, the change in W(T) 
between the points D and C is equal to i times the change in the stream function between 
these points which in turn is equal to the net flow through the nozzle. Hence, if Vj is 
the average velocity at the nozzle exit and A is the nozzle width 

iVjA = iiM-K, 0) - ii/)(K, 0) = W(-K + iO) - W(K + iO) 

Using the relations (ref. 7) 


cn(-K) = cnK = 0 


am 0 = 0 
am(2K) = n 
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in equation (9a) gives 


W(-K) = - — C 0 
k' 0 


W(K) = (tt + C ) 
k' 0 


so that 


iV T A = iA? 

1 k’ 

Hence the constant A is related to the flow by 

V T k'A 

A = -1 (10) 

7T 

As shown in sketch (d), the impingement on the plate causes the jet to divide; a por- 
tion of the net flow through the nozzle, designated by 6^, flowing to the right and a 
portion 6^ to the left. The flow to the right must be equal to the jump in the imaginary 
part of W at the point E. In order to evaluate this jump notice that 

9 m W(-K + irj) = 9m W(-K) = -Ac 

k’ 0 

lim 9mW(-K + iT]) = - — C_ (11) 

Vr> ° 

77— K* K 

On the other hand since 


cn(£ + iK’) = - L^Ll 
k sn £ 


sn(| + iK’) = 

k sn i 


tn(| + iK’) = -A- 
dn ij 
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it follows from equation (9b) that 


W(| + iK’) = - ^{sn | tan -1 (- + tan"Y— V “ + c 

k* ] 41 \ k* sn£/ \dn U 2 ° 

- — ^-isn iJk’sn I - dn d + ln [dn f_ + k' 

2k’ 1 n |_ k ' sn I + dn L dn I - k* 


- i(ir + 2C Q ) 




Since k' 2 sn 2 £ - dn 2 £ = -cn 2 £ we get 




W(£ + iK') = <sn In 
2k’ 1 11 


(k'sn £ - dn 


lk' 2 sn 2 £ - dn 2 


+ 4*Li±*q - i 
Ldn | - k'J 


i (ir + 2C Q )f 

J 


= — < sn In (k'sn £ - dn£) 2 + ln( dn * + k ' ) - sn ln(-cn 2 $) - i(ir + 2CJ 

2k’ I rt \dn g - k'/ w 


= — < 
2k' 


sn £ H ln(k’ sn | - dn £) 2 - sn | w In cn 2 £ 




m /dn { + k’\ . i(sn ^ 2C ) 

\dn 4 - k’ / a ° 


It is shown in reference 9 that dn | is strictly greater than k’ (i. e. , never equal 
to k’) for -K < £ < K. 

Also, keeping in mind that | is real, the imaginary part of W is 


9m W(| + iK’) = - (1 + sn £„) 

2k’ n k’ 


AC, 


Hence, 


AC, 


lim 9m W(| + iK') = - AlL (l + S n | R ) - — £ 


4— K 


2k' 


k' 


( 12 ) 


Since the velocity at the point E must be V , it follows from equations (11) and (12) that 
if 6 R is the asymptotic height of the righthand free streamline above the ground 
sketch (d). 
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6r v 0 = lim 
rj-K’ 


9tn W(-K + i rj) - lim 9*tt W(£ + iK') 
i— k 



» AC 

+ — (1 + sn ^tt) + 

2k’ a k’ 


= — (1 + sn 
2k’ n 


Inserting A from equation (10) 


— = — — (1 + sn 4 W ) (13) 

A V 0 2 H 


As shown in sketch (d), 6 L is the width of the stream resulting from the flow to the left. 
Continuity requirements dictate that 


6 R V o + 6 L V o = V I A 


Hence, using equation (13) 


— » — - (1 - sn ? H ) (14) 

A V 0 2 


and 


or 


6r = (1 + sn g R ) 
6j^ (1 - sn ijjj) 


R 


- 1 


sn l H - 


6 R 

— + 1 
6 t 


(15) 
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Determination of Flow Quantities 


f 


Integration to yield coordinates in physical plane. - Knowing the complex conjugate 
velocity £ and the complex potential W as functions of the parametric variable T, the 
physical variable z can be found as a function of T by using equation (1). Carrying out 
the integration in the T -plane, instead of in the W-plane, 


rr 

z = i dWdT 

Jt G +iK' S dT 


(16) 


where the origin of the physical coordinate system has been chosen at the point G. After 
substituting equations (4), (8) and (10), equation (16) becomes 


z 

A 



(k sn sn T - l)dT 
dn T f2(T; | R , k) ^/f2(T; £q, k) 


(17) 


Average velocity at nozzle exit. - In order to find an expression for Vj/V Q note 
that the nozzle lies between T = -K + iO and T = K + iO, and therefore 

z(-K) - z(K) = Ae ior 


Using equation (17) 


VI^I 
7T \ V, 


i 


' K (k sn £ H sn T - l)dT 


dn T fi(T; | R , k) ^(T; 4 G , k) 


= 1 


The integration can be performed along the line DC. Therefore 


€ 


(1 - k sn £ r sn |)d£ 
dn £ n(* ; ^ H , k) l G , k) 


-1 


(18) 
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Pressure coefficient . - The pressure coexlicient Cp is defined as 
[p(x, 0) - p q ]/1/2 pV^. By use of Bernoulli’s equation 

2 

(19) 

where the velocity ratio and corresponding physical position along the ground are given 
by equation (B13). 


C P =1 - 


u(l, K’) 


Summary of Analytical Relations 

For convenience the results of the ANALYSIS are now collected in one place 


^-1 

6 t 


sn l H = 


J R 


+ 1 


(15) 




= + H 


G T S H 


n(T; k) = 


f— (T - Q 

UC [\ 

\4K r 

2K/ 


S 4 (— (T + £ - 2K) 

4 \4K r 2K/ 


IS1 = e i[(7r/2)-a] n(T; ^ k ) y n ( T; ^ k ) 


W(T) 

V o A 


= sn £ R tan -1 ^- cn 


+ am(T + K) + C 


z(T) _ k' [ 11 ] e i a 
A 7T \V, 


f 


(k sn sn T - 1) dT 
dn T f2(T; | R , k) ^(T; | Q ,k) 


(5) 


(3) 


(4) 


(9a, 10) 


(17) 
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i min 





5 r V T 1 

— =— - (1 + sn | H ) 


A V 2 
o 


(13) 


V T 


f 

-K 


(1 - k sn (ijjSn |)d£ 


dn | £2(1; k) ^£2(4j £q> k) 


-1 


(18) 


C P =1 - 


u(4, K’) 


(19) 


These expressions are further developed into working formulas in appendix B 0 


RESULTS AND DISCUSSION 

The results that will be presented are the free streamlines, wall pressure coeffi- 
cient, velocity distribution at the nozzle exit, and the ratio of the flow amounts to the 
right and left as a consequence of the jet impingement. These quantities will be pre- 
sented in terms of three physical parameters, the tilt angle a of the nozzle, the 
horizontal distance of the nozzle center from the plane of symmetry between the two jets, 
and the vertical height of the center of the nozzle exit plane above the ground. The coor- 
dinates to the nozzle center divided by the nozzle width A, are designated by X N and Y N - 

The angle a appears explicitly in the analysis and hence could be directly specified 
as an input variable in the computer program used to evaluate the anlytical expressions. 
The nozzle position however has to be determined from the calculations and the following 
procedure was used. A value was chosen for 6j J 6 R and for k, each being between 0 
and 1. Then can be found from equation (15). For a given a, is found from 
equation (5), the value of K having been found as the complete elliptic integral, equa- 
tion (2). The average velocity from the nozzle is computed from equation (18) with the 
^ in the fl functions found from equation (Al). The width of the flow stream to the 
right 6 r /A is then computed from equation (13), and 6 l /A is found by dividing Sr/A 
by the specified ratio 6jV 6 R . 

The contours of the free streamlines were obtained from the working formulas (B14) 
and (B15) developed in appendix B. This gives the free streamline coordinates relative 
to the points C and D in sketch (d). Since S L is known, it is added to the horizontal dis- 
tance between C and F on the left free streamline to horizontally locate point C. The ver- 
tical position of C is found from 6 R , the streamline height between E and D, and the 
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nozzle tilt angle. In this manner the nozzle position is found in terms of the variables 
a, k which are the convenient input quantities in the computer program. 

The relation between 6 jV 6pj, k, and the dimensionless coordinates Xj^, Yjg (the 
location of the center of the nozzle exit plane) for various nozzle tilt angles is shown in 
figures 1(a) to (e). These plots are helpful for performing computations of the flow 
configurations in addition to those that will be presented here since they provide the 
input parameters 6jV 6^ and k for a desired nozzle position. The main utility of 
these figures is that the ratio 6jV <5j^ can be readily found for a given nozzle position 
and tilt angle. From 6-^/ 6-^ the fraction of the flow that travels back up the y-axis is 
found as 


6 L + 6 R 




+ 1 


In figure 1(a), which is for 0° nozzle tilt, as the spacing X N between nozzles is in- 
creased, the ratio SjVSr approaches 1. For large Xj^ the two jets have little influ- 
ence on each other and the flow from each jet divides equally. When the nozzles are 
tilted, for a fixed the 6 L / 6^ approaches a constant value less than 1 as the spacing 
X N is increased. The limiting values of 6jV 5^ for large X^ are shown on the 
curves; for example, for a = 15° and Yjq- = 1. 5, 6jV = 0. 588 for all Xj^ > 2. 

The velocity distribution along the nozzle exit plane, and the velocity along the 
ground were computed from equations (B12) and (B13). The velocity variation along the 
ground was inserted into equation (19) to obtain the pressure coefficient. The free 
streamlines, velocity distribution across the nozzle exit plane, and the pressure coeffi- 
cient along the ground are plotted in figures 2 to 6. 

Each of figures 2 to 6 is for a fixed tilt angle; 0°, 5°, 10°, 15°, or 20°. The (a) 
and (b) parts of each figure are respectively for spacings X N = 1 and 2, and each part 
has the two vertical positions Y N = 1 and 2. 

An examination of the figures reveals a number of specific characteristics of the 
flow. When the nozzle is near the ground, the reaction against the ground causes the 
velocity distribution across the nozzle to be nonuniform and decreases the flow out of the 
nozzle. Raising and/or tilting the nozzle decreases the interaction with the ground and 
causes the velocity distribution to become more uniform. The nozzle exit velocity is 
also influenced to a small extent by the spacing between the nozzles; the velocity dis- 
tribution being more uniform as X^. is increased. 
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The effect of nozzle tilt on the amount of flow back along the y-axis is shown in 
figure 7. Figures 7(a) and (b), respectively, show the results for 4 nozzle tilt angles at 
the fixed nozzle position X N = 1, Y N = 2, and 3 tilt angles for X N = 2, Y N = 2. In fig- 
ure 7(b) the back flow decreases with increasing tilt angle. This is the behavior that 
would be expected from intuition and does occur for most nozzle positions. Figure 7(a) 
for a smaller X^ shows the effect of a more complex interaction of the two jets. In 
this instance as the tilt increases, the back flow increases until a = 23°; then a further 
increase in a causes the back flow to decrease. This is shown more explitly in fig- 
ure 8, where 6-^/5^ has been cross plotted as a function of tilt angle. 


CONCLUDING REMARKS 

Conformal mapping was used to obtain the configuration of flow from two slot nozzles 
tilted with respect to the ground. The condition imposed is that the flow leaves each 
nozzle in a direction perpendicular to the nozzle exit plane. The flow strikes the ground 
and a portion flows back upward along the plane of symmetry between the nozzles. The 
interaction of the flow from each nozzle with the ground and with the flow from the other 
nozzle causes the velocity distribution leaving the nozzles to the nonuniform. 

The free streamlines are given for nozzle tilt angles up to 20°; for VTOL applica- 
tion only small tilt angles would be of interest. At most nozzle positions, the back flow 
between the nozzles decreased with increasing tilt angle. For small spacings between 
the nozzles an interesting effect was found. Tilting the nozzles increased the back flow 
until a maximum was reached; then with larger tilting the back flow began to decrease. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, October 30, 1968, 

129-01-07-07-22. 
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APPENDIX A 


PROPERTIES OF FUNCTION Q 


The properties of the function will be derived in this appendix. To accomplish 
this it will be useful to list some properties of the Theta -functions. The Jacobi notation 
as modified by Tannery and Molk for the Theta -functions will be used here. This is the 
same notation as used in reference 10. There are four types of Theta -functions denoted 
by Dj, $2> a 3> The functions 0^, and ^3 are defined in terms of 0^. For any 
complex constant r such that 9m t > 0, the Junction $ 4 is defined by 


Qjz |t) = 1 + 



cos 2nz 


(Al) 


for every complex number z = x + iy. It is shown in reference 10 that the series (Al) is 
absolutely convergent in the whole complex plane and that 0^ is an entire function 
(holomorphic in the whole complex plane). It is also shown in references 10 and 11 that 
the zeros of 0 4 are all simple. Thus, if z q is a zero of $ 4 

0 4 (z|t) ~ (z - z Q ) 

for z near z Q 

(A2) 

The zeros of 0 4 are the points z m n given by (ref. 11) 


= — 7tt + mir + mr 
m, n 2 

; m, n = 0, ±1, ±2, . . . 

(A3) 

It is seen from the behavior of cos 2nz 

in equation (Al) that 


0 4 (z±tt|t) 

= 3 4 (z 1 r) 

(A4) 

3 4 (-z|-r) = 

"P 

N 

ll 

(A 5) 


The remaining Theta -functions are now defined in terms of 0 4 by (ref. 10) 


0j(z| t) 


- ie lz+ (l/4)« r f, 4 



(A6) 
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,9 2 (z | t) = + | tt | (A7) 

0 3 (z|t) = 3^z + | v\tJ (A8) 

By substituting the definition (Al) into the definitions (A6) to (A8) these Theta- 
functions are given by 

OO 

«j(z I t) - 2 ]T (.i)n e l»T[M.(l/2)] 2 sln(2n + 1)2 (A9) 

n=0 

OO 

5 2 (z[t) = 2 ^ e irr[n+(l/2)] cos(2n + l)z (A10) 

n=0 


^(z |t) = 1 


+ 2 ^e 

n=l 


mm 


cos 2nz 


(All) 


From the behavior of the sine function together with eq. (A9) it can be concluded that 


^(-z |t) = -^(zlr) 


(A 12) 


For the problem analyzed in this report the argument r in the Theta -functions is 
always a pure imaginary number, and it will be assumed that r is pure imaginary in the 
discussion that follows. In this instance their definitions show that the Theta -functions 
are all real when z is real; that is, along the real axis of the z -plane. It follows from 
(A3) that there are no zeros of on the real axis, and hence it follows that 3^ cannot 
change sign there. From (Al) at z = x + iO = tt/2 

= l + 2^e i7rn2T >0 

n=l 



Since cannot change sign on the real axis 
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3 4 ( x |t) > 0 


- °° < X < + 


OO 


(A13) 


To examine the complex conjugate of 3^ and 3 3 note that cos z = cos x cosh y - 
i sin x sinh y, which yields 


(cos z)* = cos z* 


Then from equations (Al) and (All) 


[~3 4 (z J r) J *= 3 4 (z*|t) "" 

[a 3 (z|r)]* = ^ 3 (z*|r) 

Now let r be the closed rectangle in the complex T -plane shown in sketch (i) 


(A14) 


-K + iKV 


-K + iO 



K + iK' 


K + iO 


(i) Region r in T-plane showing boundary point r 
used in derivation. 


with vertices at the points -K + iK', K + iK’, K and -K. For any real number £ such 
that -K < £ r < K, fi is the complex function defined on r by 


n(T; £ r , k) 


U- (T - * r ) 

iK*\ 

\4K r 

2K / 


d 


<5 


(T + 4 r - 2K) 


Ter 


(Al 5) 


which is equation (3) in the main text. 

Since ^ is an entire function and has simple isolated zeros, it is clear that O is 
meromorphic in the interior of r (i. e. , it can have only isolated poles in F which would 
occur when in the denominator of (A15) is zero). In fact the formula (A15) can be 
used to extend O to a meromorphic function in the entire T-plane. The function S2 can 
have only simple poles located at the points T for which 
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(A16) 


3 4 (— (T + l r - 2K > 

4 \4K r 



= 0 


Equation (A3) shows that these points are given by 

T = -4 r + 2(1 + 2m)K + i(l + 2n)K’ ; m, n = 0, ±1, ±2, . . . 

However, for all Ter (and recalling that | l r | < |K|) it must be true that |£ + £ r | <2K, 
and so equation (A16) shows that has no poles in T or on the boundary of T. Thus, 

Q, is holomorphic in the interior of T, and since the poles are isolated, £2 can be 
extended to a holomorphic function in some open set which contains I\ This shows that 
S2 is continuous on the boundary of I*. 

The zeros of £2 are located at the points T for which 


(- 

\4K 


(T - 4 r ) 



= 0 


Equation (A3) shows that these points are given by 

T = l r + 4mK + i(l + 2n)K’ ; m, n = 0, ±1, ±2, . . . (A17a) 

Since |£ -4 r | < 2K and 0<tj<iK’ for TeT, the only zero of J2 on T is at a 
boundary point on the upper side of T given by (m = n = 0) 

T = £ r + iK’ (A 17b) 

and this must be a simple zero of f2. Thus, from equation (A2) 

«(T; * r ,k) ~ (T - l r - iK’) for T near l r + iK’ 

Now the behavior of S2 will be examined on the sides of the rectangle r. Equation 
(A 13) shows that on the bottom side of r which lies on the real axis of the T -plane, 

O(T;£ r ,k)>0 T€DC (A 18) 


Along the right side of T it follows by use of equation (A5) that 
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n(i ?7 + K; 4 r , k) = 




(K - |_ + il7) 


iK'\ 


2K/_ 


1 3 4 (— (K - 4 r +ir?) 

,4K 


(— (K - ! r -i T?) 
\4K r 


JjL 

2K/ \4K 


(K - + i??)* E] 



(A19) 


Then by applying equation (A 14) 


J2(i7] + K; l r , k) 


3 4 ( — (K - | + ir?) 

4 \4K 


iK'\ 
2K / 


^ 4 f— (K - 4 + ir?) 
4 \4K r 


iK’ 

2K 


This shows that along the right side CF of the rectangle T 

| J2(T; 4 r , k)| = 1 for T e CF 


(A20) 


Now examine on the left side of T where T = i7j - K. By applying equations (A4) 
and (A 5) 


Q,(ir] - K; 4 r , k) = 


* 4 | 

— (4 + K - ir?) 
y4K r 

iK'\ 
2K / 

s 4l 

f-iL (* +K + iij)* 

Uk r 

iK’\ 
2K / 

*4 

f— (4 r + K + i7?) 
\4K r 

iK* ^ 
2K , 

) 8 4l 

{— (? r + K + ir?) 

Uk r 

iK’ \ 
2K / 


(A21) 


Then by applying equation (A14) 


n(i V - K; 4 r , k) 


(4 + K + ^ 

L \4K r 


iK* 

2K 




i(-' 

\4K 


*4[^<«r +K + i *> 


iK' \ 
2K / 


Hence, 


|«(T; | , k) j = 1 for T e DE 


(A 2 2) 


On the top side of T, T = 4 + iK'. From equations (A15) and (A6) it follows that 
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(A23) 


n(S + iK’; 4 r , 


k) = 




/—(I 

\4K 


r 2 2K 


(« + 4 - 2K) +1^1 
4K r 2 2K 



iU/2)[(4 r /K)-l] 

e 



Since (— )(4 - 4 ) and + 4„ - 2K) are real, we see that 3^— (4 - f ) 

\4K/ r \4K/ r X \4K r 

5 l(— ^ + 4 r ~ 2K) 

\4K r 


iK' 

2K 


and 


iK T \ 

] are real. Hence the argument of £2(4 + iK'; 4 , k) is either 

2K / 


equal to or can differ from the argument in the exponent, — 


2 \ K 


1 , only by tt depending 


on whether or not the ratios of the i9j's is a positive or negative number. Equation 
(A17a) shows that in traversing the point | + iK’ along the contour FrE in a counter- 
clockwise direction (as shown in sketch (i) and illustrated at the points 4^ + iK’ and 
4 r + in sketch (h)) the argument of £2 must decrease by u. Since 4 r + iK' is the 
only zero of £2 on r, this is the only point of FrE where the argument of £2 can 
change. Hence, the argument of £2 is constant on Fr and on rE. Letting 4 = K in 
(A23) gives 


£2(K+ iK'; 4 r , k) = e 


_ i0r/2)[(4 r /K)-l] < K - 'M 


4K 


(K - O 

iK'\ 

X 

2K / 

(4 r - K) 

iiu\ 


2K / 


Using equation (A12) this becomes 


£2(K+ iK'; 4 r , k) = 


i(7r/2)[(4/K)-l] 

e 


i(ff/2)[(4_/K) + 1] 
e 
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Thus, 



;T £ Fr 


;T£rE 


(A24) 


In summary then we have shown that J2 is holomorphic on the interior of T and 
continuing on its boundaries. It has modulus of 1 on the vertical sides of r, and is 
positive on the bottom side of r. The argument of J2(T; | , k) on the top side of r is 
equal to (7r/2)[(£ r /K) + l] for £ > £ r and (n-/2) [(^ r /K) - l] for £ < 
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APPENDIX B 


WORKING FORMULAS FOR COMPUTING FLOW PATTERN 

In this appendix equations (4) and (17) will be rewritten in various forms which are 
convenient for computing the velocities on, and the positions of, the various boundaries of 
the flow. 

It will be convenient to make the following definitions 


£ + = — (K - 0 , 
4K 


4 = JL(K-U, T? + = — 17, 
r 4K r 4K 



R^ 4) d) = 
= 

iQ^d) 

e 

e 


s 4 U + -Ijlr 4 ) 


s 4 « + * 


• + ) 

^d + 

-41- 

r + ) 

^d + 

+ 4l' 

T + ) 

_ V4 

+ i7j + 

K 4 ) 


- i7? + 

!- + ) 

_ V4 

+ irj + 

It 4 ) 


- iV + 

K) 


(Bl) 

(B2) 

(B3) 

(B4) 

(B5) 


Equation (A13) shows that R^d) is real and positive, and equation (A9) shows that 
R^d) is real. Equations (A14) show that the ratios on the right side of equations (B4) 
and (B5) each have a modulus of 1; hence, Q^(?j) and Q^(?|) are real. 

The function O defined in equation (3) can then be evaluated along the sides of T, 
corresponding to the boundaries of the flow, by use of the functions R and Q. It follows 
from the definition of S2 and equation (A 5) that 

£>d; £ r , k) = R^ 4) d) (B6) 
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Equations (A12) and (A23) show that 


fi(|+ iK’; | r> k) »e 


ife r /K)-l]U/2) (1) 


R^'(l) 


(B7) 


Equation (A 19) shows that 


S2(K + iij; | r , k) = e r 
Equations (A21), (A 5), and (A 8 ) show that 


fi(-K + itj; £ r , k) = e 




(B 8 ) 


(B9) 


To put equation (17) in a form suitable for numerical computation purposes, it is 
necessary to express (k sn sn T - l)/dn T in terms of real quantities when T is 
complex. The following relations are deduced in reference 9: 

sn(± K + iri) = ± cn(i 7 /)/ dn(i 7 j) 

dn(± K + iri) = k'/dn(i 7 j) 

dntir], k) = dn(i], k’)/cn (rj, k') 

cn(i7j, k) = l/cn(7j, k') 


It follows that 


sn(± K + it), k) = ±l/dn ( 77 , k’) 

dn(± K + irj, k) = k' cn(t], k')/dn(? 7 , k') 

k sn(l H , k)sn(± K + irj, k) - 1 ^ dnfa, k’) + k sn(£ H , k) 

dn(± K + it], k) k' cn ( 77 , k’) 


(BIO) 


Since (ref. 9) 


sn(£ + iK’) = 1/k sn £ 
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dn(£ + iK') = -i cn £/sn | 


it follows that 


k sn sn(| + iK’) - 1 
dn(£ + iK’) 


(sn £ - sn | R ) 
cn £ 


(Bll) 


Using equations (B6) to (Bll) and equation (5) in equations (4) and (17) shows after the 
obvious deformations of the path of integrations: 

Pos itions and corresponding velocity along nozzle exit: 


z(-K) - z(g) = 
A 


'vA* C l i 
J-K Hi 


k sn(£ H , k)snd, k) 


■K dn(£ 


, k )R§ } d)^^d)~ 


d| 


> -K < 4 < K 




(B12) 


Position and corres ponding velocity along the ground : 

x({. K’)_t v l\k< I' sntt, k) - sntt H , k) 




^ / 


d G end, k)Rg^) 

Rg>(« ^>tt) 




■<d 


d) 


>-K < 4 < l G (B13) 


Location of free streamlines: 


x(-K,t?) -x(-K, 0) _ If II 
A 


1®)= ifc) f H».*> C osri- a+ Q( 3 )(,).iQ(| ( „)] d ; 

■n Y o/Jq end, k') >-2 2 J 


y(- K, rj) - y(-K, 0) _ 
A 


1 77 [ dn(T? ’ k,)+k3n(l H >k) ] sin [E 


end, k’) 


— - ar+Qjj^T?) + — Qq^(t?) 1 dr? 

L2 H 2 U J 


J 

(B14) 
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[dn(7], k’) - k sn(l H , k)] 
cn(7j, k') 

[dn(7?, k’) - k sn(£ H , k)] 
cn(77, k') 


COS 

- - a + Q$(v) + - Qq ' Hv)\&y 

_2 2 _J 



> 

sin 

i. a+Q g) W+ lQg) w j 

dij 

-J 


where 0 < rj < K' 


The following quantities are useful for the computation of 



and 


(B15) 


^ 2 

aJ. 4 ^ 4 ) = 1 + 2 E (- l) n q n cos(2n|t)cosh(2n77 + ) 
r n=l 


B^(7? + ) = -2 2^ (- l) n q n sin(2n£j)sinh(2ntj + ) 
n=l 

00 2 

A^(?? + ) = l+ 2^]q n cos(2n|^)cosh(2n7j + ) 
n=l 


y (Bi6) 


B^(?7 + ) = -2 q 11 sin(2n|^)sinh(2n7} + ) 

n=l 


where q = e 


(tt/2)(K'/K) 


Using the relation 


cos(x + iy) = cos x cosh y - i sin x sinh y 

substituting equations (Al) and (All) in definitions (B4) and (B5), respectively, and 
taking real parts yields 
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(B17) 


must always be taken in the first and second quadrants and must always be 

taken in the third and fourth quadrants. 
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(a) Coordinates of center of nozzle exit plane; X^ « 1 ; Y N = 1 and 2. 



0 .40 .80 1.20 1.60 2.00 2.40 2.80 3.20 3.60 4.00 4.40 4.80 

Dimensionless coordinate, x/A 


(b) Coordinates of center of nozzle exit plane: X N = 2; Y N = 1 and 2. 

Figure 6. - Flow configuration and wall pressure coefficient for 20° nozzle tilt angle. 
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